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Abstract 

Dirac, Schwinger and Zwanziger theories of electric and magnetic charges are obtained 
via duality transformation. Analogous construction for three Euclidean dimensions, with 
magnetic charges interacting with electric currents, is also done. The role of Dirac strings 
as dislocations in the configurations of gauge potential is emphasized. 
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1 Introduction 

In this paper we obtain Dirac, Schwinger and Zwanziger theories [1, 2] of electric and magnetic 
charges via duality transformation. Our aim is to build unified techniques for handling a field 
and its dual on an equal footing. The reason is that the dual field plays an important role in 
many contexts. It behaves as a disorder parameter and drives the properties of the theory in 
some phases. The dual field couples locally to certain topological configurations of the original 
field. Therefore by keeping them we will be able to naturally handle some non-perturbative 
aspects of the theory. A plasma or condensate of such topological configurations may be 
qualitatively and quantitatively relevant in some phases. By having a formalism which has 
both fields, we can handle such effects on the original field. Also combinations of the field and 
its dual close together often have exotic properties and play crucial role for the properties of 
the theory. 
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We illustrated our techniques for a scalar field in two Euclidean dimensions in [3]. The end 
point of a line discontinuity is the source for the dual field, much like the end point of a Dirac 
string behaving as a magnetic monopole in electrodynamics [1]. The two-dimensional local 
Lagrangian involving both the scalar field and its dual [4] is the analogue of Zwanziger's [2] 
two-potential local theory of electric and magnetic charges [5, 6]. The dislocation line becomes 
invisible for a quantization of dual charges as in Dirac theory of magnetic monopoles. The 
correlation of the field with its dual has unusual properties. 

All these show that the two-potential formalism is not restricted to electrodynamics. Any 
field theory can be recast as a local theory of the fields and their duals present together. This is 
important for quantum chromodynamics because confinement property is expected to be driven 
by topological configurations such as monopoles and vortices. We need to be able to see their 
effects on gluons and quarks. Therefore a formalism with both fields together is very useful. 

The scheme of this paper is as follows. We use Euclidean formalism throughout this paper 
to highlight the role played by >/— T- In Sec. 2, we start with a real massless scalar field in 
three Euclidean dimensions. We relate it to Abelian gauge theory by a duality transformation. 
Point sources for the scalar field are mapped to Dirac strings acting as dislocation lines in 
the configurations of the gauge potential. We further relate this to a local theory with both 
the scalar and gauge potential present simultaneously and coupling locally to the magnetic 
charges and electric currents. We demonstrate how the interactions amongst magnetic charges 
and electric currents is recovered. In Sec. 3, we begin with the Abelian gauge theory in three 
Euclidean dimensions and recover the same formalism. Here the electric currents act as sources 
of surface dislocations in the configuration of magnetic scalar potential. In Sec. 4, we obtain 
the Dirac, Schwinger and Zwanziger formulations of electric and magnetic charges via duality 
transformation in four Euclidean dimensions. We discuss the relevance of our techniques in 
Sec. 5. In the Appendix, we give some techniques useful for handling the Dirac potential of a 
magnetic monopole. 



2 Line dislocation and a local action with both the 'pho- 
ton' and the 'dual photon' in three Euclidean dimen- 



We begin with a free massless real scalar field x( x ) i n three Euclidean dimensions. Its correlation 
functions can be obtained from the functional integral 



Here JVi is a normalization factor such that Z[p — 0] = 1, and p(x) is an external source 
coupling locally to x( x )- We have deliberately included y~ T in this source term for convenience 
in performing the duality transformation below. We linearize the dependence on x( x ) i n t ne 
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exponent in equation (1) by introducing an auxiliary field B(x): 

Z[p] = N 2 J ^5(x)^x(^)e /d3:E ["^(-) 2 -^(-)-VxW+^)xW] _ ( 2 ) 

A formal integration over x( x ) gives 

Z[p]=N 3 J &B\\5{V ■B{x)+p{x))e^ ix ^ 3(x)2 ] . (3) 

X 

For a point source 

p(x) = g8 3 (x-y), (4) 

the ^-functional constraint in (3) corresponds to a magnetic monopole of strength g at the point 
y. For solving it, we choose the particular integral in the form of a Dirac string from y along 
the — ve z-direction. This singular solution gives a net flux g through any surface enclosing y. 
Thus the solution for a general p(x) is 

B(x) = V x A(x) - ii 3 d 3 1 p{x) (5) 

where 

POD 

d 3 1 P( x ) = ~ dx 3 p(x 1 ,x 2 ,x 3 ) (6) 

J X3 

and ri 3 is the unit vector in the 3-direction. The use of equation (5) in equation (3) rewrites the 
massless scalar theory of (1) as an Abelian gauge theory. The gauge field has one transverse 
degree of freedom in three Euclidean dimensions, matching that of the scalar theory. 

We shall refer to A(x) as the 'photon' and x{ x ) as the 'dual photon'. We are interested in 
their mutual correlations. We therefore include a source j(x) for A(x): 

Z[pJ]=N 4 J @Ae fd3x [~^ ^A^)-nzd^ P {x)) 2 +ij(x)-A(x)\^ ^ 

The dual photon xi x ) couples locally to the magnetic monopole density p(x). Thus equation 
(7) gives a (gauge) theory with both electric current and magnetic charges. Equation (7) shows 
that a point magnetic charge at y has the effect of a line dislocation (the Dirac string) starting 
at y in the configuration space of the gauge potential. 

As a consequence of these singular dislocation lines, the configurations A(x) which matter 
in the functional integral (7) are not the usual plane waves. For the action to be finite, V x A 
should also be singular and cancel the Dirac string singularities. Thus the configurations which 
matter are precisely the Dirac potential A D (x — y) due to a magnetic monopole at y and its 
distortions. This is explicitly seen as follows. Let us shift A(x) to a(x) as 

A(x) = a(x) + Jd 3 yA D (x-y)p(y). (8) 
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Now use the representation of A D (x — y) in the form of Green function for 93V 2 , as given in 
(54). This leads to 

V x A(x) = V x a(x) + -L / d 3 y-^-p(y) + n 3 d^p(x). (9) 

An J \x — y\ A 

So this shift cancels the Dirac string in (7). The second term on the r.h.s of (9) is simply the 
magnetic field at x due to a magnetic monopole density p{y). It is the gradient — Vx(x) of a 
scalar potential 

X (x) = J d 3 yA(x-y)p(y), (10) 

A(f-y) = ^t^- (11) 
47r \x — y\ 

Therefore 

Z[p j] — N4 J Qjae^ ' d3x [-\{^^ x )-^^ x )) 2+i 3i x )- n i x )]+ i ! d ^ ^2) 

Now 

(V x a(x) - Vx(x)) 2 = (V x a{x)) 2 - 2V • (x(x)V x a{x)) + V • { X V X ) + x{x)p(x) (13) 

as V 2 x{%) — ~p{ x )- There is no boundary contribution from the total divergence terms. So 
we get 

Z\p j] = N 4 J ^J/ d M~H^ x "W) 2+ ^^ (14) 

which has a conventional action for the new fluctuations a(x). This completes our contention 
that the configurations that contribute to (7) are the Dirac potential j d 3 yA D (x — y)p(y) and 
its fluctuations. We may fix a gauge for a(x) as usual and integrate over a(x) in (14). This 
gives 

Z[p,j] = N 5 e/^^M - !^ 1 )^ 5- ^^ ? (15) 

where A(x — y) and Ay(f- y) are respectively the propagator for a massless real scalar and 
the Abelian gauge potential in three Euclidean dimensions [7]. It shows the electric currents 
j(x) interacting via the Biot-Savart law. (Gauge fixing permits us to extend the law to currents 
that need not be conserved.) It also has the magnetic monopoles interacting via the Coulomb 
potential. In addition it shows that the magnetic monopoles p(x) interact with electric current 
j(x) through the Dirac potential A D (x — y). ( Note the ^f^l in this term in equation (15), even 
in our Euclidean theory. It is not strange as the interaction of the current j(x) with a gauge 
potential A(x) is given by ij(x) ■ A(x) even in the Euclidean theory.) This represents the net 
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effect of the line discontinuity (the Dirac string) in the configurations of the gauge potential. 
Only the end point of the string matters and behaves like a magnetic monopole. 

The magnetic charge density p(x) has a non-local coupling to the 'photon' field A(x) in 
(7), though it couples locally to the dual photon xi x ) i n (!)• We now present a local action 
that has both the photon and the dual photon fields present together. For this we rewrite (7) 
introducing an auxiliary field b(x): 

Z[p,j] = N 6 J ^&^A e ^ d M _ ^(z) 2 +jb(z)-(vx (16) 

Thus p(x) couples locally to d^b^x), which is to be identified with the dual photon x( x ) : 

X (x) = 9s%(x) . (17) 
We may integrate back over bi(x) and b2(x), to get 

Z\p J} = N 7 J $ x $Ae$ dAx \-~' 2{d ' ix{x))2 ~^^ (18) 

This gives the local field theory of electric currents and magnetic charges in three Euclidean 
dimensions. It is the analogue of the two-potential formalism in four dimensions [2] and of 
the local field theory involving the scalar field and its dual in two dimensions [3, 4]. Note the 
following unusual features: 

• The action is not manifestly rotation invariant. Nevertheless, the rotation covariance is 
restored for physical observables when the Dirac quantization condition for electric and 
magnetic charges is met. (See below.) 

• The 'kinetic energy' terms for x( x ) an d A{x) have derivatives only in some directions. 

— * 

However, if we integrate over x( x ) (correspondingly A(x)), we recover the conventional 
action for A (x) (correspondingly x( x ))- 

• The action in (18) is not real (even with the imaginary sources switched off). The term 
bilinear in x an d A(x) is pure imaginary. 

The 'propagators' can be calculated using Fourier modes, as follows. In the action of 
equation (18), we use the identity (n 3 x (V x A)) 2 = (V x A) 2 - (n 3 • (V x A)) 2 . The 
resulting — |(V x A) 2 term in the action becomes \A ■ V 2 A ) upon adding the gauge-fixing term 

— * —* — * — * — * — * 

corresponding to the Feynman gauge. Also n 3 • (V x A) = (n 3 x V) • A, which is ik± ■ A{k) in 

— * — * 

momentum space, where k± = n 3 x k. Thus we obtain 

Z[pJ] = N 7 J @x(k)@A(k) exp J d 3 k( - l - x (-k)h 2 X (k) - 1 -A l (~k)(k 2 5 lJ - k ±l k ±] )A 3 (k) 
+iA i (-k)k u k 3X (k) + iji(-k)Ai(k) + ip(-k) X (kj) . (19) 
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The propagators are then obtained by the inversion of a matrix: 



—ik?k i 



In position space, the propagators are 



1 



k 3 



° k 3 



(x(x)x(y)) 



A(x-y) 

Aij(x-y) = 6ijA(x-y) 
iA°{x-y) 



(20) 



(21) 
(22) 
(23) 



where A(x — y) is given by (11), and we have used (57). Alternatively, we can read off these 
position space propagators from (15). 

Even though x( x ) an d A(x) are real fields, the propagator (Ai(x)x(y)) is pure imaginary. 



This is possible because the action is not real. We see that %A 



d, 



x 



y) serves as the 'propa- 



gator' connecting the electric currents and magnetic charges. The correlation of x( x ) with the 

— » — * — * 

'magnetic field' B(x) = V x A(x) has the Dirac string singularity: 



(B(x) X (0)) 



1 x 

47T \x\ 3 



+ n 3 S(x 1 )S(x 2 )0(-x 3 ] 



(24) 



Because of the explicit presence of the Dirac string, rotation invariance in (23) and (15) is not 
manifest. Dirac [1] argued that with a 'quantization' of electric (e) and magnetic (g) charges, the 
Dirac string becomes invisible and rotation covariance is restored. Consider a point magnetic 
charge given by (4), and a loop C carrying a current 



dT dXdr) 5 ^ 
dr 



At)) 



(25) 



where r is an arbitrary parametrization of the loop C. The contribution to the cross-correlation 
of exp \ie § dx l A i ( y x)\ with exp [igxiu)} comes from the last term in the exponent of (15) for the 
sources (4) and (25). (This is the analogue of the cross-correlation between the vertex operators 
for the scalar field and the dual field in two dimensions, which was considered in [3].) Thus 
this cross-correlation equals exp [ieg § dx 1 Af \x — y)\ . Using Stokes' theorem, this is 



exp 



-i^-n(c) 

47T 



(26) 



where fl(C) is the solid angle subtended by C at the site y of the magnetic charge, and the 
solid angle is to be computed by using a surface bounding C which does not intersect the 
Dirac string. Therefore when an infinitesimal loop C does not enclose the string, we get the 
contribution 1, but when it encloses the string, we get the contribution e~ ieg . Only with the 
quantization condition 

eg = 2im (27) 
the latter contribution is also 1, and the Dirac string is then invisible to any current loop. 
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3 Surface dislocations and scalar potential theory of elec- 
tric currents and magnetic charges in three Euclidean 
dimensions 

In Sec. 2, we started with a massless real scalar field and obtained the two-potential theory 
of magnetic charges interacting with electric currents in three Euclidean dimensions. In this 
Section, we begin with Abelian gauge theory and obtain the same two-potential formalism. 
This exercise is instructive for the case of four Euclidean dimensions. 

We begin with 

Z[j] = N 8 J ^A{x)e^ Ax H^ xA{x))2+i ' {x) - A ^ (28) 

describing the current j(x) interacting via the gauge potential A(x) in three Euclidean dimen- 
sions. Rewriting 

Z[j] = Ng J ^^A e / d3x ["I^) 2 +^)-^ x ^) + ^)-^( :E )] (29) 

= N w J ®b\\ 5(V x b(x) + ](x))e^ ix H i{x)2 ] . (30) 

x 

The consistency of the constraint requires 

V-j(x) = 0. (31) 

Choosing the solution 

H x ) = diX{x) - esud^jiix) (32) 
for the 5-functional constraint, we get 

Z[p,j] = N n J ^xe /d3x [~^ (9ix(x) ~^ lj2(x))2 ~^^ (33) 

where we now introduced the source for x- Note that the component js(x) is not explicitly 
present. However, by the conservation law (31), we can write 

j 3 (x) = -d^jiix) + d 2 j 2 {x)) (34) 

and therefore it is implicitly present. 

Equation (33) is giving the interaction of magnetic charges and electric currents using the 
scalar potential xi x ) encountered in magnetostatics. Consider a current loop C in the 1-2 plane 
with a charge e flowing in it. Equation (33) shows that for the action to be finite in this case, 
there should be a discontinuity in the scalar potential x i n the form of a surface dislocation. 
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This dislocation is along a cylindrical domain wall with C as the mouth and extending all the 
way to infinity in the 3-direction. The gradient of the potential jumps by e across the domain 
wall. This is the conventional description of using a multivalued magnetic scalar potential in 
the presence of electric currents [8]. 

We linearize (33) in a specific way: 
Z[pJ] = N 12 J ® X ®A{x)\[^Az{x) - a 3 (x)) exp^ [ - ^(d 3X (x)) 2 - l -{d 2 A 3 {x) - d 3 A 2 (x)) 

-^{daA^x) - d 1 A 3 (x)) 2 - i(d 2 A 3 (x) - d 3 A 2 (x))(d lX (x) - d^j 2 (x)) 

-i(d 3 A 1 (x)-d 1 A 3 (x))(d 2X (x)+d^j 1 (x))+ip(x) X (x)\(35) 

Equation (33) is reproduced from equation (35) by shifting Ai, I = 1, 2 to A'j = Aj — d 3 ~ 1 diA 3 
and integrating over A\. The potential A 3 has been introduced in this step to have a gauge- 
invariant Lagrangian. The gauge-fixing condition A 3 [x) = a 3 (x) corresponding to the axial 
gauge is here the simplest choice as it is not affected by the shift from Aj to A\. 

Finally it can be checked, using (34), that the action in (35) is equal to the action in (18). 
The Lagrangian being gauge-invariant, we can now pass from the axial gauge to any other 
gauge we find convenient. 

Instead of choosing an infinite domain wall, we can simply choose a finite surface S enclosing 
the current loop C to be the dislocation for the scalar potential. For that case, in the place of 
(32), we have 

bi(x) = dix{x) + diXs{x), (36) 
where Xs in discontinuous across the surface S: 

V X s(x)= [ d 2 X(s)n(s)e5 3 (x-X(s)), (37) 

X(s) being a point on the surface S, and n(s) the normal to S at this point. Because of this 
discontinuity, (<9j<9j — djdi)x s (x) ^ 0, and (36) gives V x b(x) +j(x) = with j(x) given by (25). 
(This can be checked by using the identity J dS x Vip = § dltp where if) is a scalar function.) 
When (36) is used in (30), the finiteness of the action requires x to have a discontinuity across 
S so as to cancel the discontinuity in Xs- Therefore for a loop C linked to C, we get 

e ig$ c ,dx l dix(x) _ e iegN cc , ^gg-j 

where N C c' is the linking number: the number of times loop C winds around loop C in the 
clockwise sense. (This is because the integral of the current j over any open surface bounding 
the loop C equals eN C c-) Thus, (38) equals 1 when the quantization condition (27) is satisfied. 
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4 Two-potential theory of electric and magnetic charges 
in four Euclidean dimensions 



Consider the quantized Abelian gauge field in four Euclidean dimensions: 

Z\jfj\ = N 13 J S$A^e$ d4x \~^ 9 ^ A v( x )- d » A ^ x )) 2+i in A ^ x )\ . 

Here (/i — 1, 2, 3, 4) is the external current. We have 



(39) 



Z[ h ] = N u I 9b i 9e i 9A ll ex.Y> / d 4 x - ^{xf - ^{xf + iei{x)(d A A t (x) - d^x)) 



+ ieijkbiixjdjA^x) + iji(x)Ai(x) + ij A (x)A A (x) 
N 15 / @bi@ei Y[ S(diei(x) + j A (x)) JJ 5{e ijk djb k (x) - d^x) + ji{x)) 

X X 

x exp / d A x 



1 1 

- ^(x) 2 - -bi(x) 2 



From the 5-functional constraints we have the consistency condition 

diji(x) + d A j A (x) = 0. 
We solve the first constraint in (41) as 

ei(x) = e ijk djC k {x) - S^d^j^x) , 



(40) 



(41) 



(42) 



(43) 



corresponding to choosing the Dirac string along the 3-direction (similar to (5)). Using this, 
the second constraint becomes 



eijkdj(b k - d A C k ) = -ji - 5 i3 d 3 1 d A j A . 
We solve this in the form (similar to (32)) 

b k (x) - d A C k (x) = -d k C A {x) - erndz 1 ji{x) 

Therefore we get 



(44) 



(45) 



= N 16 / ^C M exp / d 



x 



-(diC 2 (x)-d 2 C 1 (x)-d 3 1 j 4 (x)f 



-(djCsix) - d 3 Cj(x)) 2 - \{diCi{x) - d A d(x) + eud^jjix)) 2 
\{d 3 C A {x) - d A C 3 (x)) 2 + tk^x)C,(x)] 



(46) 
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where we have introduced a source k^ for C M . Here the indices /, J run over only 1 and 2. This 
gives the Dirac and Schwinger formulations. As in Sec. 4, we linearize (46) in a particular way. 
The first term in the Lagrangian is linearized to 

- l -{d 3 A 4 {x) - d 4 A 3 {x)f + i(d 3 A 4 (x) - 9 4 A 3 (x))(9 1 C 2 (x) - d 2 C 1 (x) - d 3 l ]4 (x)) 
while the third term is linearized to 

-^(d 3 Mx) - djA 3 {x)) 2 + ieuidsAjix) - djA 3 {x)){djC 4 {x) - d 4 Cj(x) + e JK d 3 l j K {x)) . 

The axial gauge condition is imposed on A 3 . (We get back (46) by shifting A 4 and Af, and 
then integrating.) This linearized form of (46), on using (42), is equal to 

Zlj^k,] = N 17 J ^C^A.HdiA^-a^exp J d 4 x[-^(d I A 3 (x)-d 3 A I (x)) 2 

X 

-\{d 3 A 4 {x) - d 4 A 3 (x)) 2 - ^(djC 3 (x) - d 3 Cj(x)) 2 - l -{d 3 C 4 {x) - d 4 C 3 (x)) 2 
+t(d 3 A A (x) - 9 4 A 3 (x))(9 1 C 2 (x) - d 2 C 1 (x)) 

+ie IJ {d 3 A I {x) - djA^idjC^x) - d 4 Cj{x)) + ij^x)A^{x) + ik^C^x) 

(47) 

Thus we have recovered Zwanziger's two-potential theory of electric and magnetic changes via 
a duality transformation. If we had chosen h as the direction of the Dirac string instead of the 
— ve z-axis, we would have got for the exponent on the r.h.s of (47) 



d 4 x 



\{n ■ E{x)f -Unx B(x)f - \{h ■ £{x)f - \{h x B(x)) 2 + i(n ■ B{x)){n ■ E{x)) 



+i(h x S(x)) ■ (n x B(x)) + ij^A^x) + ik^(x)C^(x) 



(48) 



where B = V x A, E = -VA 4 + d 4 A, B = V x C and S = -VC 4 + d 4 C. 



5 Discussion 

For a variety of reasons, it is useful to have a formulation with both a field and its dual field 
simultaneously present in a local theory. In [3] this was done for a scalar theory in two Euclidean 
dimensions and the advantages were highlighted. In this paper, we have carried this out for 
Abelian gauge theory in three and four Euclidean dimensions. 

The general features are: 

• The sources for a field are certain types of singular dislocations in the configurations of 
the dual field and also vice versa. 
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• The role of these dislocations is to force discontinuous boundary conditions on the fields. 
Thereby new sectors of the field configurations are explored. 

• A local theory with both the field and its dual present simultaneously has certain unusual 
features. Though there are more fields, it is equivalent to the original theory and the 
degrees of freedom are not changed. This happens because the dual fields are hidden 
in the auxiliary fields as specific non-local combinations. As a consequence correlations 
of fields with their duals have unusual properties. The theory is not manifestly rotation 
invariant. However, rotation covariance is recovered for the 'right' observables with a 
quantization of the charges of the field and the dual field. 

• These features are already known in the context of Dirac's theory of magnetic monopoles. 
Our thrust is that they are general properties of dual fields and not restricted to elec- 
trodynamics. We have obtained Dirac, Schwinger and Zwanziger formulations of electric 
and magnetic charges via duality transformations. We have emphasised the role of Dirac 
string as singular dislocation in the configurations of the gauge potentials. These issues 
are relevant for non-Abelian gauge theory. Many non-perturbative aspects such as con- 
finement are expected to be driven by topological configurations which couple locally to 
the dual field. This will be discussed elsewhere. 
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A Appendix 

In this Appendix, we represent the Dirac vector potential of a monopole in the form of a Green 
function. The Dirac potential of a monopole located at the origin has the form 

?r>, s 1 sm # -? 1 - ? /^x 

A (x) = -0 = — n 3 x 49 

v ' 4nr(l + cos6>) 4tt r + x 3 K ' 

with the Dirac string along the — ve z-direction. Here r = \x\ and f = xj\x\.. [For checking 
(49) and other results below, a useful formula is n 3 = cos Of — sin 99.] Let us write 

A D (x) = n 3 x c (50) 

where the vector field c is undetermined upto addition of a vector in the 3-direction. We choose 

c = ~a (51) 

Att r + x 3 
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so that [9] 



c = V/, / = ^ln(r + x 3 ). (52) 



Note that d 3 f = l/47rr, and so d 3 V 2 f = — Thus the Dirac potential at x due to a 
monopole at x' can be expressed in terms of the Green function for the operator <9 3 V 2 : 

A D (x — x') = — n 3 x V ln(|a; — x'\ + x% — x' 3 ) (53) 

4:71 

= -n 3 x VPaVW-x')]. (54) 

[An alternative form of the Dirac potential is A D (x) = — 0(l/47rr) cot = — h 3 x f (x 3 /47rp 2 ), 
with the Dirac strings along the ±z directions. Here p 2 = x 2 + x 2 2 . This alternative form 
is half of the sum of the two forms of A D (x), one having the Dirac string in the — ve z- 
direction and the other having the Dirac string in the +ve z-direction. In this case, we choose 
c— (rh 3 — x^/Anp 2 in Eq. (50). Then c — V/ and 83 f = l/Anr continue to hold, but with 
/ = (I/871") ln((r + xs)/(r — X3)). So equation (54) is still valid] 

The result given in Eq. (54) can also be seen by going over to the Fourier space. For the 
potential of Eq. (49), 

V x A D (x) = ^ + n 3 5(x 1 )5(x 2 )9(-x 3 ) . (55) 
Taking the Fourier transform, we get 

kxA°$) = -± + *±. (56) 
(The Fourier transform of the theta function can be obtained using d9(x)/dx = S(x).) We now 

— * — > — * 

evaluate kx both sides and use k ■ A D (k) = (since V • A u (x) = 0) to obtain 



AO C k) = |^ . (57) 



This agrees with Eq. (54). 
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